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ABSTRACT 

Tracking the spectral evolution of transiently accreting neutron stars between outburst and quies¬ 
cence probes relatively poorly understood accretion regimes. Such studies are challenging because they 
require frequent monitoring of sources with luminosities below the thresholds of current all-sky X-ray 
monitors. We present the analysis of over 30 observations of the neutron star low-mass X-ray binary 
SAX J1750.8-2900 taken across four years with the X-ray telescope aboard Swift. We find spectral 
softening with decreasing luminosity both on long (~1 year) and short (~days to week) timescales. 

As the luminosity decreases from 4 x 10 36 erg s _1 to ~ 1 x 10 35 erg s _1 (0.5-10 keV), the power law 
photon index increases from from 1.4 to 2.9. Although not statistically required, our spectral fits allow 
an additional soft component that displays a decreasing temperature as the luminosity decreases from 
4 x 10 36 to 6 x 10 34 erg s _1 . Spectral softening exhibited by SAX J1750.8-2900 is consistent both 
with accretion emission whose spectral shape steepens with decreasing luminosity and also with being 
dominated by a changing soft component, possibly associated with accretion onto the neutron star 
surface, as the luminosity declines. 


1. INTRODUCTION 

Low-mass X-ray binaries (LMXBs) are powerful and 
strongly variable X-ray sources that provide a means to study 
accretion flows over a wide luminosity range and on rela¬ 
tively short timescales. Neutron star LMXBs consist of a neu¬ 
tron star accreting via Roche-lobe overflow from its low-mass 
(< 1M 0 ) companion. In transient LMXBs, the accretion rate 
onto the neutron star surface varies by many orders of mag¬ 
nitude, but the systems spend the majority of their time in 
a low-luminosity, 10 31 — 10 33 erg s~ 4 (0.5-10 keV unabsorbed 
luminosity, denoted by Lx), state of quiescence during which 
there is little or no accretion. In outburst, these systems un¬ 
dergo intense accretion for periods of weeks to months (to 
years in the extreme case of quasi-persistent transients), with 
typical X-ray luminosities of 10 36 —10 38 erg s _1 . In this work, 
we focus on the intermediate luminosity range between out¬ 
burst and quiescence of 10 34 — 10 36 erg s _1 . 

Limited instrumental sensitivity and/or observation time 
makes it difficult to study neutron star transients during 
their, often rapid (lasting less than several weeks), tran¬ 
sition from outburst to quiescence. Hence, the accretion 
processes and emission sources in this luminosity regime 
are not well-understood. However, as more and more neu¬ 
tron star binaries are being monitored between outburst 
and quiescence, there is growing evidence that these sys¬ 
tems may share the common behavior of spectral softening 
with decreasing luminosity. Spectral softening between out¬ 
burst and quiesence has been documented in detail for tran¬ 
sient black hole binaries (Yuan et al. 2008; Wu & Gu 2008; 
Plotkin et al. 2013) and for individual neutron star transients 
(Armas Padilla et al. 2011, 2013; Linares et al. 2014). It has 
recently been proposed as a feature exhibited by most neu¬ 
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tron star transients between luminosities of 10 36 and 10 34 erg 
s' 1 (Wijnands et al. 2014). In black hole systems, the soften¬ 
ing is commonly attributed to radiatively inefficient accretion 
flows. The underlying mechanism(s) responsible for the soft¬ 
ening in neutron star systems is likely more complicated due 
to the presence of the neutron star surface and magnetic field. 

A better understanding of accretion at luminosities between 
10 34 3 4 and 10 36 erg s _1 may provide insight into accretion pro¬ 
cesses at even lower luminosities while the systems are in 
quiescence. In addition to thermal emission from the neu¬ 
tron star surface, many quiescent transients display a non- 
thermal, high-energy spectral tail above 2 keV which can be 
modeled with a power law (r = 1-2) 5 . Nonthermal and ther¬ 
mal emission may be variable both in their total flux and rel¬ 
ative flux contribution over short and long timescales, as has 
been found in the well-studied neutron star transient Cen X-4 
(Cackett et al. 2010). Several quiescent transients have also 
exhibited flares, sudden increases in flux reaching ten or more 
times the typical quiescent emission levels and exponentially 
decaying over the course of a few days (Fridriksson et al. 2011; 
Degenaar & Wijnands 2013). The physical origin of the non¬ 
thermal emission in quiescence is not entirely understood, but 
is generally attributed to ongoing accretion (Campana et al. 
1998), while flares are interpreted as sporadic increases in the 
accretion rate onto the neutron star, although the physics of 
disk instabilities at low accretion rates is not well-understood. 

More frequent and deeper observations of neutron star 
LMXBs between their quiescent and outburst states, as well 
as observations of flares in quiescence, are critical to con¬ 
straining the physical origins of the hard and soft spectral 
contributions, together with their evolution in time. Such 
work is necessary to understand the nature of accretion at 
all luminosities. In particular, it can help in identifying the 
accretion regimes that heat the neutron star crust, which has 
important implications for equilibrium temperature measure¬ 
ments of the crust and studies of the thermal evolution of 
transient neutron stars. 

SAX J1750.8-2900 was first detected in outburst in 1997 
by the Wide Field Cameras on board BeppoSAX and has 
since been detected in outburst an additional 3 times with 

5 where T is defined as dN/dE oc E~ r 
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both long (4-5 months) and short (< 1 month) outburst du¬ 
rations. The presence of type I X-ray bursts revealed that the 
source is a neutron star LMXB (Natalucci et al. 1999). Using 
RXTE observations of its 2001 outburst, the second recorded 
for this source, Galloway et al. (2008) analyzed photospheric 
radius expansion bursts to establish a distance of 6.79 ± 0.14 
kpc for hydrogen-poor burning. In March 2008, SAX J1750.8- 
2900 entered its third known outburst (Markwardt & Swank 
2008) with a reported flux of 230 mCrab corresponding to a 
luminosity 3.1 x 10 3 ' erg s _1 (2-10 keV, d = 6.79 kpc), simi¬ 
lar to the peak flux exhibited in previous outbursts. By Au¬ 
gust 2008 the source fell below the RXTE detection threshold 
and was reported to be returning to quiescence (Linares et al. 
2008a). Super AGILE reported a 60 seconds long burst on 
October 9, 2008 (Pacciani et al. 2008) indicating the system 
had rebrightened. The return to outburst was also seen in two 
Swift XRT observations on October 1 st and 8 th (Linares et al. 
2008b). By mid-February 2009, SAX J1750.8-2900 had re¬ 
turned to quiescence (Lowell et al. 2012). 

Over a year after the end of its last reported outburst, 
SAX J1750.8-2900 was observed for the first time in quies¬ 
cence in April 2010 with ATMMwith a bolometric luminosity 
of 1.05 ± 0.12 x 10 34 (D/ 6.79 kpc) 2 erg s' 1 ( L x = 8.9 x 10 33 
erg s _1 ), making it the most luminous quiescent neutron star 
known (Lowell et al. 2012). The spectrum was well-fit with 
a neutron star atmosphere model with a high effective tem¬ 
perature (kT^g = 148 eV). Under the assumption the source 
had undergone no additional heating (i.e. accretion) since 
outburst, Lowell et al. (2012) claimed that the neutron star 
crust and core were in equilibrium and, thus, the high surface 
temperature reflected a hot neutron star core. 

In the two years after the XMM observation, SAX J1750.8- 
2900 exhibited a short outburst (< 1 month) and a flare in 
quiescence. A February 2011 increase in activity reported by 
IBIS/INTEGRAL (Fiocchi et al. 2011) was confirmed as a 
faint outburst by Swift XRT with Lx = 1-05 x 10 36 erg s^ 1 
for a source distance of 6.79 kpc (Natalucci et al. 2011). In 
March 2012, the detection by Swift XRT of a low-level flare 
(with characteristic Lx = (3 — 4) x 10 34 erg s _1 assuming 
d = 6.79 kpc) lasting between 5 and 16 days indicated the 
system may undergo small accretion events between outbursts 
(Wijnands & Degenaar 2013). 

No orbital parameters are known for this neutron star bi¬ 
nary and no observations of the companion star exist except 
for reports of a near-IR candidate counterpart (Torres et al. 
2008). The line-of-sight towards SAX J1750.8-2900 suffers 
from high-extinction (Ay 6 = 13.0-15.1), consistent with the 
high column density, Nh > 10 22 , which may limit observa¬ 
tions of the optical counterpart. 

In this work we present a detailed spectral analysis of the 
neutron star LMXB SAX J1750.8-2900, observed with Swift 
XRT over the course of four years. We have tracked the X- 
ray spectral evolution of SAX J1750.8-2900 over three orders 
of magnitude in luminosity between outburst and quiescence, 
detecting spectral softening with decreasing luminosity. We 
interpret the spectral softening as being generated by changes 
in soft and hard accretion-powered emission as the mass ac¬ 
cretion rate decreases, and we investigate the potential crustal 
heating effects due to accretion at sub-outburst rates. 

2. OBSERVATIONS & ANALYSIS 

We analyzed all Swift XRT (Burrows et al. 2005) observa¬ 
tions where SAX J1750.8-2900 was in the field of view (FOV), 
which yielded 27 source detections out of a total of 38 obser¬ 
vations. Five observations had SAX J1750.8-2900 more than 
7 arcminutes off-axis (i.e., inside but near the edge of the 

6 http://irsa.ipac.caltech.edu/applications/DUST/ 
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Figure 1. Light curves for SAX J1750.8-2900 from RXTE PCA 
(top) and the Swift BAT monitor (bottom). In April 2010, indi¬ 
cated by the diamond symbol, SAX J1750.8-2900 was observed in 
quiescence with XMM and Lowell et al. (2012) reported a quies¬ 
cent luminosity of Lx = 9 x 10 33 erg s — L The red vertical line 
symbols in the bottom panel indicate all 27 Swift XRT observa¬ 
tions with source detections, analyzed herein. Within this set of 
observations, SAX J1750.8-2900 had luminosities between 10 34 and 
10 37 erg s _1 (see Table 2). 

FOV), but we included them in our analysis as luminosity 
upper limits. The data set spans over four years of activity, 
between March 2008 and September 2012, with a total expo¬ 
sure time of almost 70 ks. The long-term soft and hard X-ray 
light curves from SAX J1750.8-2900 measured with RXTE 
PCA and Swift BAT (Krimm et al. 2013), respectively, cov¬ 
ering the timespan of our Swift observations are shown in 
Figure 1. The majority of our Swift XRT observations cover 
the tail end of the 2008 outburst and the return to outburst 
just several weeks later in October; the remaining observa¬ 
tions primarily occurred in February and March 2012. A log 
of the Swift XRT observations is presented in Table 1. 

2.1. Source Detection & Spectral Extraction 

All observations were re-processed with the xrtpipeline 
tool (v.0.12.6). The source was identified in photon counting 
(PC) mode observations with the sosta tool set to a signal- 
to-noise threshold limit of 2.0 in Xlmage (v.4.5.1). 

Pile-up corrections were made for PC mode observations 
with count rates greater than 0.5 s _1 within a circular region 
with a 20-pixel (47 arcsec) radius centered on the source. 
The extent of the piled-up region was determined following 
the Swift XRT Pile-Up thread'. In these observations, an 
annulus with an inner radius set to exclude the piled-up pixels 
(radius of 3 pixels for mild pile-up and up to 12 pixels for the 
most severe case of pile-up, ObsID 31166001) and an outer 
radius of 30 pixels (45 pixels for 31166001) was used for the 
source region. The background spectra were extracted from 
two circular regions, both with radii of 45 pixels. All PC 
mode observations with count rates between 0.004 and 0.5 
s^ 1 (i.e., not affected by pile-up) had spectra extracted from 
a circular region with a 20-pixel radius centered on the source, 
and the background was extracted from two circular regions, 
both with radii of 35 pixels. Potentially nearby X-ray sources 
(within 15 arcseconds) were identified using Simbad 8 and were 
avoided in any of the source and background regions. 

Window timing (WT) mode observations only incur pile- 

7 http://www.swift.ac.uk/analysis/xrt/pileup.php 

8 http: //simbad.u-strasbg.fr/simbad/ 
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Table 1 

Exposures and background corrected count rates for all 
observations with and without SAX J1750.8-2900 detections. 95% 
upper count rate limits (indicated by the ’<’ symbol) are given in 
cases where the source was not detected above the SNR threshold 
( < 2). The indicates the far off-axis observations (> 7’) with 
SAX J1750.8-2900 within the FOV. 


ObsID 

Observation 

Date 

Mode 

Exposure 

(ks) 

Count Rate 
(s- 1 ) 

31166001 

2008-03-17 

PC 

1.0 

3.36 ±0.06 

31174001 

2008-03-24 

PC 

0.7 

1.56 ±0.05 

31174001 

2008-03-24 

WT 

0.3 

5.71 ± 0.14 

31174002 

2008-03-25 

PC 

0.9 

1.38 ± 0.04 

31174002 

2008-03-25 

WT 

0.1 

3.96 ± 0.23 

31174003 

2008-08-14 

PC 

1.9 

0.06 ± 0.01 

31174004 

2008-08-15 

PC 

1.9 

0.04 ±0.01 

31174005 

2008-08-16 

PC 

2.0 

0.06 ± 0.01 

31174006 

2008-08-17 

PC 

2.1 

0.05 ± 0.01 

31174007 

2008-08-18 

PC 

2.1 

0.10 ±0.01 

31174008 

2008-08-19 

PC 

1.7 

0.12 ± 0.01 

31174009 

2008-08-20 

PC 

2.0 

0.08 ± 0.01 

31174010 

2008-08-24 

PC 

1.8 

0.06 ± 0.01 

31174011 

2008-08-25 

PC 

2.1 

0.13 ± 0.01 

31174012 

2008-08-26 

PC 

2.5 

0.12 ± 0.01 

31174013 

2008-08-27 

PC 

1.9 

0.11 ± 0.01 

31174014 

2008-08-28 

PC 

1.7 

0.08 ±0.01 

31174015 

2008-08-29 

PC 

5.9 

0.050 ± 0.0003 

31174016 

2008-09-03 

PC 

2.3 

0.14 ±0.01 

31174017 

2008-09-05 

PC 

0.7 

0.05 ± 0.01 

31174018 

2008-09-11 

PC 

1.8 

0.07 ±0.01 

31174019 

2008-10-01 

PC 

2.0 

0.74 ± 0.02 

31174020 

2008-10-08 

PC 

1.9 

0.62 ± 0.02 

31174021 

2011-02-19 

PC 

3.9 

0.46 ± 0.01 

31174023 

2011-09-30 

PC 

0.5 

<0.018 

31174024 

2012-02-14 

PC 

3.8 

<0.004 

31174025 

2012-02-26 

PC 

2.6 

<0.003 

31174026 

2012-02-29 

PC 

0.3 

<0.023 

31174027 

2012-03-03 

PC 

3.2 

0.002 ± 0.001 

31174028 

2012-03-06 

PC 

2.8 

<0.003 

31174030 

2012-03-17 

PC 

3.1 

0.018 ±0.003 

31174031 

2012-03-20 

PC 

1.0 

0.006 ± 0.003 

31174032 

2012-03-22 

PC 

1.0 

<0.006 

43649001 

2012-05-22 

PC 

0.5 

<0.019* 

32427001 

2012-06-17 

PC 

0.8 

<0.009* 

32427002 

2012-06-22 

PC 

0.7 

<0.008* 

32427003 

2012-08-15 

PC 

0.1 

<0.042* 

32427004 

2012-09-04 

PC 

0.5 

<0.010* 

Quiescence 

- 

PC 

12.3 

0.0011 ± 0.0004 

Flare 

- 

PC 

5.1 

0.012 ± 0.002 


up for count rates above 100 s^ 1 ; all SAX J1750.8-2900 count 
rates were below 10 s^ 1 . Thus, no pile-up corrections were 
required in this mode. Both the source and background spec¬ 
tra in WT mode were extracted from circular regions with 
30-pixel radii. 

After extracting the source and background spectra in XS- 
elect (v.2.4b) with standard event grades 0-2 for WT and 
0-12 for PC modes, we created exposure maps and ancillary 
files using the xrtexpomap and xrtmkarf tools, respectively, 
using the most recent RMFs (v.14 for PC mode and v.15 for 
WT mode) and enabling the vignetting and PSF correction 
options. Whenever possible, we grouped energy bins to a min¬ 
imum of 20 counts per bin using the grppha tool. In cases 
where there were only dozens of counts, we binned as low as 
5 counts per bin. We confirmed that the spectral parameters 
obtained with binning to fewer than 20 counts per bin were 
consistent with the spectral parameters from fits of unbinned 
spectra. 
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Figure 2. SAX J1750.8-2900’s count rate light curve: The down¬ 
ward arrows indicate the upper 95% count rate limits for all ob¬ 
servations without source detections (see Section 2.1.1), and the 
inset shows the activity between August and October 2008 where 
the bulk of our observations occur. The numbers indicate the ap¬ 
proximate groups of observations with similar spectral properties 
and luminosities. 
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2.1.1. Non-detections and combined observations 

In 11 out of a total of 38 observations, we did not detect our 
source. We were able to place upper limits on SAX J1750.8- 
2900’s activity using a circular extraction region with a radius 
of 15 pixels centered on the published Chandra coordinates 
for SAX J1750.8-2900 (Chakrabarty et al. 2008). Following 
the prescription of Gehrels (1986), we placed the 95% upper 
limits on the source count rate (denoted by < in Table 1). 

We combined observations with no source detection and ob¬ 
servations with faint source detections at low intensities (with 
count rates < 0.02 s -1 ) during the February and March 2012 
period to increase the source counts and perform spectral 
fits. Observations 31174030, 31174031 and 31174032, cor¬ 
responding to March 17-22, exhibited elevated count rates 
compared to previously exhibited quiescent levels (< 0.002 
s -1 ); we combined these three observations to form a ’Flare’ 
spectrum. We formed a ’Quiescent’ spectrum by summing all 
observations immediately prior to the flare with count rates 
(including upper limits) less than 0.004 s -1 for a total ex¬ 
posure time of 12.3 ks (corresponding to ObsIDs 31174024- 
25 and 31174027-28). Spectra for the flare and quiescence 
were generated by summing observations’ event files in XS- 
elect followed by summing their exposure maps in Xlmage. 
The source and background spectra were extracted following 
the same procedure previously outlined for PC mode observa¬ 
tions with clear source detections using a circle with a 10-pixel 
radius as the source region. 

The source intensity as a function of the observation date 
for our entire data set is shown in Figure 2. 

2.2. Spectral Analysis 

All spectral analysis was performed using XSpec v. 12.7.0 
(Arnaud 1996). We used the tbabs multiplicative model to 
account for interstellar absorption and set our abundances to 
those of Wilms et al. (2000) and used the cross-sections from 
Verner et al. (1996). All quoted luminosities are the 0.5-10 
keV unabsorbed luminosities (unless a different energy range 
is specified), calculated with the cflux convolution model 
assuming a source distance of 6.79 kpc (Galloway et al. 2008). 
The Eddington fraction (Lx/L ecLcl) assumes an Eddington 
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luminosity, L^dd, of 2 x 10 38 erg s _1 for a 1.4 Mq neutron star 
with a hydrogen-rich photosphere. Errors on fit parameters 
correspond to the 90% confidence bounds (A\ 2 = 2.71) and 
all chi-squares, \ 2 ( dof ), are reduced. 

2.2.1. Hardness Ratio 

In order to study the evolution of the average spectral slope 
or ’’hardness” in a model-independent way, we calculated the 
hardness ratio (HR) for each observation. We have defined 
the hardness ratio as the ratio of hard (2-10 keV) to soft 
counts (0.5-2 keV), which we calculated within XSpec from 
background-corrected source spectra binned to 1 count per 
bin to avoid having empty bins. 

Swift XRT analysis guides 9 warn of soft spectral ’bump’ 
below 1 keV for highly absorbed sources ( Nh > 10 22 cm -2 ) 
when using the standard WT mode event grades 0-2 selection. 
We investigated the grades 0-2 versus grade 0 extractions for 
our two WT mode observations and found no significant dif¬ 
ference in the hardness ratios between the two grade selections 
(ObsID 31174001 HR = 5.7 ±0.4 vs 5.95 ±0.5, 31174002 HR 
= 6.6±1.3 vs 6.1 ± 1.1 for grades 0 and 0-2, respectively) or in 
parameters obtained from spectral fits (for example, a simple 
absorbed power law). Thus, we used the standard WT event 
grades 0-2 extraction in our analysis. 

2.2.2. Column Density Constraints and Spectral Modeling 

The column density towards SAX J1750.8-2900 is poorly 
constrained; published values range from Nh = (2.5 — 
6.0) x 10 22 cm -2 , henceforth noted as Nh, 22 = 2.5 — 6.0 
(Natalucci et al. 1999; Lowell et al. 2012) and are strongly de¬ 
pendent on the spectral model used. From a quiescent XMM 
observation, Lowell et al. (2012) obtained Nh ,22 = 5.9 ± 0.5 
with the NSATM0S model (Heinke et al. 2006), while a black- 
body fit yielded Nh, 22 = 4.0lo g. 

Fitting individual observations with an absorbed power law 
with no fixed inputs led to large uncertainties in the spec¬ 
tral parameters. In order to obtain better constraints, we 
fit groups of observations with a single model and tied cer¬ 
tain parameters between all observations within the group. 
We formed the groups by initially fitting all individual obser¬ 
vations with an absorbed power law fixing the column den¬ 
sity to Nh, 22 = 5 (the average of the values obtained by 
Lowell et al. 2012) while leaving the photon index and nor¬ 
malization free to vary; we then sorted the observations by 
their 2-10 keV unabsorbed luminosities 10 . Plotting the 2-10 
keV unabsorbed luminosities against the powerlaw photon in¬ 
dex, we found 7 distinct clusters of observations in the photon 
index-luminosity parameter space and divided our set of ob¬ 
servations into 7 spectral and luminosity groups, shown in 
Table 2. 

We also attempted a single component fit with a black- 
body model (bbodyrad), but found it did not provide accept¬ 
able fits across our dataset’s entire luminosity range. At high 
luminosities {Lx > 10 37 erg s , the blackbody yielded a 
significantly worse fit than a power law (y 2 = 1.2 vs 0.9). 
Additionally, column densities required for the blackbody fits 
{Nh, 22 < 2.5) were lower than even the low-end of published 
estimates. 

9 http://www.swift.ac.uk/analysis/xrt/digest_cal.php 

10 We first attempted to sort the observations by their 0.5-10 keV 
unabsorbcd luminosity, but found a weaker correlation between 
the luminosity and spectral shape (i.e. photon index). For highly 
absorbed sources {Nh ,22 ± 3) and soft spectra (r > 2), a power 
law model leads to large 0.5-2 keV unabsorbed fluxes due to the 
divergence of the power law flux at low energies. Thus, we found a 
stronger spectral dependence on luminosity when we excluded the 
low-energy flux contribution by looking at the 2-10 keV unabsorbed 
luminosity. 
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Figure 3. The 0.5-10 keV unabsorbed luminosity estimates for 
SAX J1750.8-2900 between February and March 2012, including 
the flare reported by Wijnands & Degenaar (2013) which begins 
as early as March 6. Both luminosity upper limits (indicated by 
the downward pointing arrows) and luminosity estimates for ob¬ 
servations with source detections (diamond symbols) were calcu¬ 
lated from the count rate-flux conversion based on a blackbody 
with = 331 eV. The shaded region is bounded by the qui¬ 

escent luminosities reported by Lowell et al. (2012) for blackbody 
{Lx = 2.8 X 10 33 erg s —1 and neutron star atmosphere models 
{Lx = 8.9 x 10 33 erg s _1 ). 

We performed a new power law fit for each group, tying 
the photon index while leaving the normalization free. For 
the two brightest groups (6 and 7) we left the absorption 
parameter free (but tied if there were multiple observations in 
the group) and found consistent values for the column density, 
Nh ,22 = 4.2 ± 0.4 and 4.4 ± 0.4. We fixed all future column 
densities to Nh ,22 = 4.3, the average value from our two 
brightest groups. 

While the group fits provided tight constraints on SAX 
J1750.8-2900’s average spectral changes with luminosity, to 
examine the source’s variability on shorter timescales we fit 
the individual observations with an absorbed power law, fix¬ 
ing the column density while leaving the photon index and 
normalization free to vary. Due to insufficient source counts, 
we were unable perform spectral fits for ObsIDs 31174024- 
31174032 in an attempt to track the changes in the photon 
index before, during and after the flare. Instead, we used a 
count rate-flux conversion to estimate the system’s luminos¬ 
ity. Using WebPimms (v4.6a) n , we converted the observed 
count rate/upper limit to a luminosity based on the pub¬ 
lished blackbody quiescent fit result of Lowell et al. (2012) 
{kTf§ = 331 eV and N H , 22 = 4), shown in Figure 3. The 
quiescent and flare spectra formed by combining observations 
(Section 2.1.1), however, had enough counts to perform power 
law fits. 

While the power law model worked well for both the group 
and individual fits (we obtained \ 2 ~ 1), we investigated 
whether multiple spectral components were present by per¬ 
forming a multicomponent model fit to the grouped obser¬ 
vations where our statistics were good enough to constrain 
additional parameters. We added a blackbody component 
(bbodyrad) to the existing powerlaw. When tying the black¬ 
body temperature and normalization across all observations 
within a group, we were able to constrain the soft compo¬ 
nent’s temperature, emitting area (i.e. the blackbody model’s 
normalization) and its fractional contribution to the total un- 

11 http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 







5 



Figure 4. Hardness ratio as a function of count rate: As the 
count rate decreases between 1.5 — 0.05 s -1 , the hardness ratio also 
decreases, indicating the source softens as its intensity declines. 
This behavior is intrinsic to the source as the hardness ratio is 
independent of a spectral model. The source may soften further 
towards quiescence (down to a count rate of 0.001 s —3 ) but we 
discuss the limitations of interpreting the softening with low count 
spectra in Section 3.1.1. 

absorbed flux for groups 3-6. 

3. RESULTS 

In our set of observations, SAX J1750.8-2900’s count rate 
and luminosity span over 3 orders of magnitude with ranges 
of 0.001 to 6 s' 1 and Lx between 10 34 to 10 37 erg s' 1 . Ob¬ 
served at the beginning of the March 2008 outburst with 
Lx = (4 — 30) x 10 36 erg s' 1 , the source declined in in¬ 
tensity and luminosity in August 2008 (Lx ~ 1 x 10 35 erg 
s' 1 ) until it re-brightened in October 2008. In March 2012, 
SAX J1750.8-2900’s activity was below the RXTE PCA and 
Swift BAT sensitivity thresholds, but its quiescent activity 
(< io 34 erg s" 1 ) was still marginally detectable with Swift 
XRT. Our observations were all well-fit by a pure power law, 
and in some cases, we were able to fit groups of observations 
with a power law plus blackbody model. 


10 


10 


l x /l e 

10' 3 


10 " 


10 4 



Figure 5. Model 1, Power Law: Hardness ratios (model- 
independent) and photon indices from a power law fit versus the 
unabsorbed 0.5-10 keV luminosity, also given in terms of the Ed¬ 
dington fraction (Lx / LEdd)- Empty symbols indicate values for 
the individual observations and the symbol color represents its lu¬ 
minosity group. Black filled circles are the results for the fits of the 
grouped observations (with the group numbers adjacent). Spectral 
softening, first evident in the hardness ratio-count rate relation in 
Figure 4, is further supported by the anti-correlation and correla¬ 
tion of the hardness ratio and photon index, respectively, with the 
luminosity between 4 X 10 36 and 10 35 erg s~ 1 . 


3.1. Spectral Softening, Hardening and Variability 

We find a clear softening towards lower luminosities of the 
X-ray spectrum of SAX J1750.8-2900. Hardness ratios and 
photon indices decrease and increase, respectively, as the sys¬ 
tem’s unabsorbed luminosity decreases. 

In Figure 4, as the source’s intensity decreases, 1.5 -A 0.05 
counts s' 1 , the spectrum softens as the hardness ratio de¬ 
creases by nearly a factor of 3, from 8 —> 3. The results 
presented in Figure 4 use hardness ratios and count rates for 
pile-up corrected spectra and are therefore independent of the 
spectral model used. Further softening at even lower count 
rates (< 0.05 s' 1 ) is discussed in the next section. 

3.1.1. Model 1: Power law behavior 

The softening behavior is also evident in the power law 
spectral parameters both when observations are fit as groups 
and individually. Power law parameters obtained from group 
fits are shown in Table 2 and plotted in Figure 5, while power 
law fits to individual spectra are contained in Table 3 and 
plotted in Figures 5 and 6 for the full data set and the August 
to October 2008 period, respectively. As seen in Figure 5, 
as the 0.5-10 keV unabsorbed luminosity decreases, Lx = 
4 x 10 36 -4 1 x 10 35 erg s' 1 , the photon index increases from 
1.4 to 2.9. 


Below an intensity of 0.05 s' 1 and a luminosity of ~ 5 x 10 34 
erg s' 1 , the spectral parameters and hardness ratios suffer 
from large uncertainties due to the low number of counts, 
making the softening trend less obvious. Lowell et al. (2012) 
and Wijnands & Degenaar (2013) found that SAX J1750.8- 
2900’s 2010 and 2012 quiescent spectra, respectively, were 
well-fit with a thermal spectral model (either NSATM0S or 
bbodyrad). In particular, for the deep XMM quiescent obser¬ 
vation where more stringent limits on the nonthermal compo¬ 
nent could be placed, no additional power law was necessary 
and a power law component had a maximum contribution of 
4% to the total flux (Lowell et al. 2012). The thermal nature 
of the quiescent spectrum combined with our low hardness 
ratio (HR = 1.1) suggest it is possible that the softening con¬ 
tinues towards quiescence near Lx ~ 10 34 erg s' 1 . The flare, 
however, appears inconsistent with the softening behavior. It 
exhibits a higher hardness ratio and lower photon index than 
the trend would predict given the event’s luminosity, but the 
error bars are large for both parameters, so we are unable 
to claim that the flare is an exception to the softening trend 
with such limited statistics 

In addition to the general trend of softening as the luminos¬ 
ity decreases over an order of magnitude, we found softening 
occurred on very short timescales (on the order of days to 
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Table 2 

Power Law and Power Law + Blackbody Fits - Grouped Observations 





Model 1: tbabs 

* powerlaw 




Model 2: tbabs * (powerlaw + bbodyrad) 




Group 

Low a 

High a 

N ff ,22 

r 

L X 

X 2 (dof) 

n H,22 

r 

kT c 

BB-Norm 

FPL* 

Fbb d 

X 2 (dof) 

F-Test e 

HR 

1 

0.001 

0.01 

4.3 

4.1 ± 1.6 

0 ll -1-0,40 

U11 _0.08 

15.1 (18) b 




- 





1.1 ± 0.8 

2 

0.01 

0.1 

4.3 

1.7 ± 0.5 

n i « + ^.05 
°' 16 Io.03 

68.4 (59) b 




- 





5.3 ± 2.1 

3 

0.1 

0.3 

4.3 

2.9 ± 0.2 

u y -0.21 

0.9 (108) 

4.3 

o k + 0.8 
2 -1.2 

0.49+ 015 

6 2+ 8 ' 1 0 6 +0 ' 2 

° -4.9 u -0.1 

66 ± 32 % 

34 ± 32 % 

0.9 (106) 

2e-3 

3.1 ± 0.8 

4 

0.3 

5 

4.3 

2.2 ± 0.1 

1 26+ 017 
0.16 

1.0 (144) 

4.3 

1 6+ 0 - 4 
x -0.4 

0 55 + 0 06 

U b -0.05 

„ _+3.5 . n +0.1 
6 ' 7 —3.1 10 -0.1 

62 ± 12 % 

38 ± 12 % 

0.8 (142) 

7e-9 

4.6 ± 0.9 

5 

5 

20 

4.3 

1.6 ± 0.1 

15 42~^ bb 
i o. tz — o.66 

1.0 (231) 

4.3 

1 6+ 0 - 2 
±b -0.2 

1 04+ 0 - 17 
±U4 —0.15 

4 4+ 3 ' 4 14 4+ 0 ' 7 

*-1.9 14 - 4 -0.7 

79 ± 7 % 

21 ± 7 % 

0.8 (229) 

5e" 2 

7.5 ± 0.6 

6 

20 

100 

4.4 ± 0.4 

1.4 ± 0.1 

44 57+ a ‘ ti ^ 

° -3.34 

1.0 (216) 

4.3 

i 0+0.3 
XD —0.3 

1 89"*" U ' /4 
-0.52 

1 7 +i.8 44 g+3.5 

A -1.5 **-3.0 

77 ± 11 % 

23 ±11 % 

1.0 (215) 

Se-6 

6.9 ± 0.8 

7 

100 

500 

4.2 ± 0.4 

1.6 ± 0.1 

9EO 1 0+19.75 
zoo.id_ 16 0Q 

0.9 (139) 




- 





6.2 ± 0.3 


a Group luminosity bounds in units 10 35 erg s 1 for the 2-10 keV unabsorbed luminosity obtained with a power law fit. 
b C-statistic 

c Blackbody temperature (keV) 

^Percent contributions of the power law and blackbody components to the total 0.5-10 keV unabsorbed luminosity. 
e F-test probability that the added component (i.e. blackbody) improves the fit by chance. 


weeks), as well as spectral hardening as the luminosity in¬ 
creased on similar timescales. While the 2-10 keV luminosity 
was more strongly variable than the 0.5-10 keV luminosity 
(see Figure 6), between August and October 2008, which cor¬ 
responds to groups 3 and 4, hardness ratios and photon in¬ 
dices were correlated and anti-correlated, respectively, with 
the 0.5-10 keV luminosity 12 which varied between 5.9 x 10 34 
and 1.6 x 10 35 erg s _1 , 

For example, SAX J1750.8-2900 softened between Au¬ 
gust 25-29 (ObsIDs 31174011-4015); as the intensity, hard¬ 
ness ratio (6.3 —> 3), and 0.5-10 keV luminosity decreased 
(Lx = (1.42 —> 0.96) x 10 35 erg s _1 ) while the photon in¬ 
dex increased (2.0 —» 2.9). Immediately after the softening, 
the spectrum hardened; between August 29 - September 3 
(ObsIDs 31174015-4016), SAX J1750.8-2900’s luminosity in¬ 
creased to 1.6 x 10 35 erg s , while the hardness ratio in¬ 
creased to 5.4 and the photon index decreased to F = 2.0. 


3.2. Model 2: Soft emission and thermal contributions 

We were able to fit a soft component (modeled with a 
blackbody) in addition to the hard emission (modeled with a 
powerlaw) when fitting several of the groups of observations. 
There was no change in the chi-squared value compared to 
the pure power law fit, indicating there is no strong model 
preference. The multicomponent fit parameters are shown in 
Table 2 and spectra with their group fits (either powerlaw or 
powerlaw plus blackbody models) are shown in Figure 7. 

Despite group 7 being the brightest observation, it had 
fewer source counts than any of groups 3, 5 or 6. We were 
unable to calculate the lower error bound on the blackbody 
temperature when we performed a powerlaw plus blackbody 
fit, indicating the blackbody component was not definitively 
detected in this observation. For a blackbody normalization 
of 5 (an average of the normalizations found in fits on groups 
3-6), the blackbody had a temperature upper limit of 2.45 
keV and maximal flux contribution of 25%. 

Luminosity-photon index relations for the power law and 
the power law plus blackbody models are compared in the top 
panel of Figure 8. With the addition of the soft component, 
between 4 x 10 36 and 1 x 10 35 erg s~ 4 the photon index was 
consistent with being constant (r ~ 1.6). At a luminosity of 
Lx = 6 x If) 34 erg s _1 , however, the power law component 
softened (T = 2.5+ l j l '2) although the error bars are large and 
still consistent with the photon index at higher luminosities. 

Shown in the middle panel of Figure 8, the blackbody com¬ 
ponent exhibited a decreasing temperature (1.9 —> 0.5 keV, 
plotted in red) at lower luminosities (Lx ~ 4 x 10 36 —> 6 x 10 34 


12 As mentioned in 2.2.2, the suppressed variation in the 0.5- 
10 keV unabsorbed luminosity may be due to the large 0.5-2 keV 
unabsorbed flux contribution in a power law fit to highly absorbed, 
soft spectra. 



190 210 230 250 


Days since MJD 54500 

Figure 6. Model 1, Power Law: Between August and Octo¬ 
ber 2008, the hardness ratios and photon indices obtained with 
a power law model show signs of correlation and anti-correlation, 
respectively, with the 0.5-10 keV unabsorbed luminosity, indicat¬ 
ing the spectrum softens (hardens) as the luminosity decreases (in¬ 
creases) on timescales from days-weeks. Also plotted in the bottom 
panel is the 2-10 keV unabsorbed luminosity, which shows stronger 
variation in time. 


erg s _1 ). The blackbody’s contribution to the total flux (i.e. 
the thermal flux fraction, plotted in black) tends towards 
higher fractions at lower luminosities, although the uncertain¬ 
ties are large. We could only fit for an additional soft com¬ 
ponent in groups 3-6, but we have also plotted the quiescent 
thermal fraction (with a lower limit of 96% as the quiescent 
spectrum was thermal with a maximum power law contribu¬ 
tion of 4%) and luminosity from the Lowell et al. (2012) re¬ 
sults. We estimated the lower limit to the thermal fraction of 
the flare by simply dividing the quiescent thermal luminosity 
by the flare luminosity. 




































7 



Energy (keV) 




Energy (keV) 



Energy (keV) 


Figure 7. Swift XRT folded spectra of SAX ,11750.8-2900 used in this work, together with our best-fit model functions. Observations 
were combined into 7 groups (Sec. 2.2, Table 2). Groups 1, 2 and 7 are fit with an absorbed power law, while with groups 3-6 we were able 
to constrain an additional soft component, and so we plot their powerlaw + bbodyrad fits. Adding a soft component improved residuals 
slightly, but overall the spectra and their residuals are similar between the single and multicomponent fits. Groups 1 and 2 are plotted 
together but are fit separately and binned for display only. 


In the bottom panel of Figure 8, we plot the power law and 
blackbody components’ 0.5-10 keV luminosities as a function 
of the total 0.5-10 keV luminosity. When we performed a 
simple linear fit, we found the power law flux declines more 
rapidly than the thermal flux (a slope of 1.1 versus 0.70 for 
the blackbody flux), which is consistent with the thermal 
flux fraction increasing towards lower luminosities. If we ex¬ 
press each component’s luminosity in terms of the quiescent 
thermal luminosity (right axis, bottom panel), we find that 
for groups 3 and 4 the thermal luminosity is 2-4 times the 
quiescent luminosity, indicating that the quiescent emission 
contributes significantly to the thermal luminosity (approxi¬ 


mately 50% and 25% of the thermal luminosity at 6 x 10 34 
erg s' 1 and 1 x 10 35 erg s , respectively). 

4. DISCUSSION 

Our set of Swift XRT observations of SAX J1750.8-2900 
span over 4 years and 3 orders of magnitude in luminosity, 
which allows us to explore a wide range of accretion rates 
and timescales. We find that the source softens between out¬ 
burst and quiescence, Lx = 10~ 2 —> 10 -3 Land the soft¬ 
ening mechanism generates spectral variability on multi-day 
timescales. 
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Table 3 

Model 1: Power Law Fits - Individual Observations 





tbabs * 

powerlaw 



Group 

ObsID 

Nff,22 

r 

T a 

LX 

X 2 (dof) 

HR 

1 

Quiescence 

4.3 

4.1 ± 1.6 

0 n+°- 4f > 

U - J - - 0.08 

15.1 (18) b 

1.1± 0.8 

2 

Flare 

4.3 

1.7 ± 0.5 

0 16 +UU5 

U-±u -0.03 

68.4 (59) b 

5.3 ± 2.1 

3 

31174004 

4.3 

2.7 ± 0.5 

0 59+ u -^ y 

u.oy_ 0 20 

1.3 (15) 

4.4 ± 1.4 


31174005 

4.3 

3.1 ± 0.5 

1 16+ 0 - 66 
' ' LU —0.41 

0.8 (21) 

2.4 ± 0.5 


31174015 

4.3 

2.9 ± 0.2 

0 96 +0 ' 23 

u.»u_o 19 

1.5 (13) 

3.0 ± 0.4 


31174006 

4.3 

3.0 ± 0.4 

0 97+ 0 - 42 
‘ -0.31 

1.2 (21) 

3.3 ± 0.8 


31174017 

4.3 

2.7 ± 1.0 

0 96’ 1 ’ 1 ' 05 

U.UO_o 42 

0.2 (5) 

2.9 ± 1.2 


31174003 

4.3 

3.0 ± 0.5 

i OO+0.67 
1 - zz -0.40 

0.6 (9) 

2.2 ± 0.5 


31174010 

4.3 

2.9 ± 0.4 

1 07"* -0 ' 49 

1 ' u ' -0.33 

0.6 (18) 

3.4 ± 0.8 

4 

31174018 

4.3 

2.4 ± 0.4 

1 02" 1-0 ' 31 

l.UZ_o 23 

0.5 (10) 

4.6 ± 1.2 


31174009 

4.3 

2.5 ± 0.3 

i -,5+0.31 
1 - 1O -0.24 

1.0 (15) 

4.1 ± 0.8 


31174014 

4.3 

2.3 ± 0.4 

0 96 +0 ' 27 
u.»U —o 20 

0.7 (11) 

3.7 ± 0.8 


31174007 

4.3 

2.5 ± 0.3 

1 43+ 0 - 35 
—0.27 

0.9 (12) 

3.6 ± 0.6 


31174013 

4.3 

2.4 ± 0.3 

i KQ+O- 34 
l.OU_ 0 27 

1.2 (12) 

4.1 ± 0.7 


31174012 

4.3 

2.3 ± 0.2 

1 42+ 0 - 24 
^•^-0.20 

1.5 (28) 

5.5 ± 0.9 


31174008 

4.3 

2.2 ± 0.3 

1 45+0.30 
J --^°-0.24 

0.4 (18) 

4.1 ± 0.7 


31174011 

4.3 

2.0 ± 0.2 

1 49+ 0 - 20 

0.9 (16) 

6.3 ± 1.2 


31174016 

4.3 

2.0 ± 0.2 

i rr+0.20 

l.OO 0 18 

1.0 (14) 

5.4 ± 0.8 

5 

31174021 

4.3 

1.5 ± 0.1 

10.20)+“ 

1.0 (80) 

7.5 ± 0.5 


31174019 

4.3 

1.5 ± 0.1 

15-79i“;£ 

1.2 (67) 

8.5 ± 0.7 


31174020 

4.3 

1.6 ± 0.1 

18.68j+»? 

0.6 (53) 

6.6 ± 0.5 

6 

31174002 PU 

4.3 

1.2 ± 0.1 

33.49++ 

0.8 (47) 

7.5 ± 0.7 


31174002 WT 

4.3 

1.3 ± 0.3 

37.51+™ 

0.9 (18) 

6.1 ± 1.1 


31174001 PC 

4.3 

1.3 ± 0.1 

48.33+!! 

1.2 (36) 

8.1 ± 0.8 


31174001 WT 

4.3 

1.4 ± 0.1 

51-33+1! 

1.0 (94) 

6.0 ± 0.5 

7 

31166001 

4 2 +u ' 4 

-0.4 

1.6 ± 0.1 253.13+jg'oo 

0.9 (139) 

6.2 ± 0.3 


a Unabsorbed 0.5-10 keV luminosity in units of 10 35 erg s 1 
b C-statistic 


4.1. Spectral softening between outburst and quiescence 

Spectral softening is a nearly universal behavior in black 
hole transient systems as they enter quiescence (Wu & Gu 
2008; Plotkin et al. 2013) and now a trend seen across dozens 
of neutron star transients (Wijnands et al. 2014). In black 
hole binaries, one interpretation of the spectral softening dur¬ 
ing the transition from the hard state to quiescence involves 
the changing opacity of an inner advection-dominated ac¬ 
cretion flow (ADAF), although jets may also play a role in 
the softening near quiescence. The anti-correlation between 
the photon index and luminosity typically spans Lx /Lgff = 
10 -2 — > 10~ J for ~ 10 Mq black holes, while below Lx/Lffff 
= 10 —5 , there is evidence the spectrum does not soften further 
into quiescence but that the photon index remains constant 
(Plotkin et al. 2013). 

While both neutron star and black hole transients ex¬ 
hibit increasing photon indices towards lower luminosities, 
the softening mechanism(s) are not necessarily expected to 
be the same, largely due to the differences in emission sources 
between the two types of transients. Gathering observa¬ 
tions of neutron star transients with luminosities between 
10 34 — 10 36 erg s^ 1 and comparing to the photon index- 
luminosity relation for black holes found by Plotkin et al. 
(2013), Wijnands et al. (2014) found that softening black hole 
and neutron star transients are distinctly separate popula¬ 
tions. Black holes are significantly harder (r < 2) than neu¬ 
tron star transients (2 < V < 3) below 10 35 erg s . 

With the exception of accretion disk emission, black hole 
transients are almost entirely sources of nonthermal emission 
and can be modeled with a pure power law at luminosities 
below an Eddington fraction of 1%. Neutron star LMXBs, 
however, have both nonthermal and thermal emission, al¬ 
though it is unclear which component dominates in different 



Figure 8. Model 2, Power law + Blackbody: The top panel 
shows the difference in the power law photon index behavior be¬ 
tween the two models: pure power law (plotted in grey) and power 
law plus blackbody (black). When a soft component (i.e. black- 
body) is included in the fit, the photon index is essentially constant 
between Lx = 4 X 10 36 —> 10 35 erg s —1 . In the middle panel are 
the thermal fraction (black) and blackbody temperature (red) as a 
function of the total 0.5-10 keV unabsorbed luminosity. The soft 
component’s temperature is lower at lower luminosities, which is 
expected if the soft component is associated with accretion onto 
the neutron star surface. Bottom panel: Luminosity from each 
of the power law (blue) and blackbody (red) components in the 
power law plus blackbody model plotted along with a linear fit, 
along with the Lowell et al. (2012) NSATM0S thermal luminosity (at 
~ 10 34 erg s —4 ). The power law flux declines more rapidly than 
the flux from the thermal component. 

luminosity and accretion regimes. Due to the neutron star’s 
solid surface, there are multiple sources of soft X-ray emis¬ 
sion across a range of luminosities: boundary layer emission 
(Lin et al. 2007), thermal surface emission due to heat gener¬ 
ated by pycnonuclear reactions deep within the neutron star 
crust (Brown et al. 1998), and a blackbody-like spectrum due 
to low-level accretion (Zampieri et al. 1995), while accretion 
can produce nonthermal emission at almost any luminosity. 
Changes in multiple or just one of these sources can contribute 
to the spectral softening. Attempts to track the thermal and 
nonthermal contributions as a function of time and luminos¬ 
ity are further complicated when there is insufficient spectral 
quality to constrain both components. Despite these obsta¬ 
cles, there is growing evidence suggesting that changes in a 
soft component’s properties is largely responsible for spectral 
softening (Armas Padilla et al. 2013; Bahramian et al. 2014; 
Wijnands et al. 2014) and that the component is powered by 
accretion rather than heat release from the neutron star crust, 
while there may be simultaneous changes in the nonthermal 
emission associated with accretion. 


Blackbody Temperature (keV) 
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Softening has been observed in neutron star LMXBs over 
luminosities of 10 36 —X 10 34 erg s _1 , similar to the range 
where we have detected softening in SAX J1750.8-2900 be¬ 
tween outburst and quiescence. In Figure 9 we have plotted 
the photon indices and luminosities for SAX J1750.8-2900’s 
power law fits, along with the published power law parame¬ 
ters for several other NS LMXBs that have exhibited spec¬ 
tral softening. For a more complete comparison of spectral 
softening in a dozen neutron star transients, see Figure 1 in 
Wijnands et al. (2014). 

4.1.1. Soft Component’s Role in Softening 

In most of the neutron star softening sources, a power law 
provided adequate, if not good, fits to most of the observa¬ 
tions. However, at lower luminosities (< 10 35 erg s^ 1 ) and 
with higher SNR spectra, many sources required an addi¬ 
tional soft spectral component, indicating the spectra are not 
entirely non-thermal, and the softening cannot be accurately 
modeled by a power law with a steepening photon index, as 
is the case with softening black hole transients. 

In deep XMM observations of neutron star transients XTE 
J1709-267 and IGR ,11794-3030, Degenaar et al. (2013) and 
Armas Padilla et al. (2013) achieved the sensitivity to con¬ 
strain both hard and soft spectral components; both sources 
exhibited a decaying thermal temperature and constant pho¬ 
ton index while the total flux decreased over the course of 
approximately 10 hours. Due to the similarities between the 
two systems’ behaviors, Armas Padilla et al. (2013) credited 
the spectral changes in both sources to a decrease in the mass 
accretion rate onto the neutron star surface which produced 
a soft component with a decreasing temperature, while the 
nonthermal flux due to accretion also decreases. 

With SAX J1750.8-2900 when adding a blackbody, we 
found almost no change in the chi-squared values, indicat¬ 
ing there was no model preference between a pure power law 
and the power law plus blackbody and that the additional 
soft component was not statistically required. Thus, the soft¬ 
ening in SAX J1750.8-2900 from our data can be accurately 
described by a steepening photon index towards lower lumi¬ 
nosities. However, it is worth noting that when we included 
a blackbody component, we find the soft component’s be¬ 
havior is similar to that seen in other softening neutron star 
transients where additional soft components are statistically 
required, and the behavior is consistent with the interpreta¬ 
tion that the soft component is generated by accretion onto 
the neutron star. In SAX J1750.8-2900, as the total luminos¬ 
ity decreased, the blackbody displayed lower temperatures 
and its flux contribution increased, indicating the strength 
of the component was increasing relative to the hard compo¬ 
nent. Between Lx ~ 10 36 and 10 35 erg s _1 , the photon index 
of the power law component was essentially constant reveal¬ 
ing that an increasingly dominant thermal component may 
contribute significantly to the softening. We are unable to 
comment on the behavior of the power law component’s pho¬ 
ton index above Lx ~ 5 x 10 36 erg s _1 and below ~ 5 x 10 34 
erg s , as we either were unable to constrain both power law 
and blackbody components or the photon index suffered from 
large uncertainties. 

4.1.2. Variability and Spectral Hardening 

Rapid, non-monotonic variability, as well as alternating 
spectral softening and hardening with decreasing and increas¬ 
ing luminosity seen in SAX J1750.8-2900 strongly support the 
theory that accretion generates the spectral softening. With 
our dense set of observations between August and October 
2008 (Figure 6), we found SAX J1750.8-2900’s luminosity and 
spectral hardness varied on timescales as short as a day; such 
rapid variability can be generated by changes in the accretion 
rate onto the neutron star. Additionally, SAX J1750.8-2900 



Figure 9. Power law photon indices versus the unabsorbed 0.5- 
10 keV luminosity obtained from pure power law fits for multi¬ 
ple neutron star LMXBs that exhibit spectral softening with de¬ 
creasing luminosity, or hardening with increasing luminosity in the 
case of Swift J174805.3-244637 (Armas Padilla et al. 2013, 2011; 
Degenaar et al. 2013; Linares et al. 2014; Bahramian et al. 2014). 

displays alternating spectral hardening and softening over the 
course of days and weeks, which has not been seen in other 
neutron stars on such timescales. 

Spectral hardening with increasing intensity has been ob¬ 
served in the neutron star LMXB Swift J174805.3-244637 as 
its rise to the hard state was monitored with Swift XRT 
(Bahramian et al. 2014). As its luminosity increased from 
L x = 4 x 10 34 -4 10 36 erg s _1 (~ 2 x 10 -4 ->4x 10 -3 
L Edd) the photon index decreased (r = 2.6 —X 1.7) when fit 
with a pure power law. Using more complex models with a 
soft component, Bahramian et al. (2014) found an increasing 
thermal temperature with increasing luminosity, which can¬ 
not be produced via a cooling neutron star crust since the 
observations are at the onset outburst. While we could not 
track the rapid day-to-day variations in the soft component 
as SAX J1750.8-2900 hardened, the similarities between SAX 
J1750.8-2900 and Swift J174805.3-244637 further support the 
interpretation that the spectral changes are produced by hard 
and soft emission powered by a varying mass accretion rate 
onto the neutron star surface. 

While we found SAX J1750.8-2900 to be variable during 
August to October 2008, its average luminosity was approxi¬ 
mately 10 35 erg s _1 , which shows that underluminous accre¬ 
tion flows (Lx/Lvdd < 0.1%) can be stable under certain con¬ 
ditions. Namely, SAX J1750.8-2900 remained at Lx ~ 10 35 
erg s^ 1 for over one month (Figure 6) in contrast with the 
rapid decay to quiescence (< 1 week) seen in other neu¬ 
tron star transients (Aql X-l, Campana et al. (2014); SAX 
J1810.8-2069 Linares et al. (2007); 4U 1608-52 Linares et al. 
(2009)). 

4.2. Accretion and the Origin of the Soft Component 

In many of the existing discussions of softening neutron star 
LMXBs, emission associated with accretion onto the neutron 
star surface has been cited as the possible soft component 
detected in the spectra. The spectrum generated by spher¬ 
ical accretion onto an unmagnetized neutron star has been 
a long-standing issue. Zel’dovich & Shakura (1969) demon¬ 
strated that freely, radially infalling ions onto a neutron star 
produced a blackbody spectrum with a high-energy tail due 
to Comptonization with luminosities of ~ 10 35 — 10 38 erg 
s _1 . Similar results were obtained by subsequent studies even 
when additional physics were included. Zampieri et al. (1995) 
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extended simulations to even lower mass accretion rates down 
to luminosities of ~ 10 31 — 10 32 erg s” 1 and found the emer¬ 
gent spectra were significantly hardened compared to a black- 
body at the neutron star effective temperature. Deufel et al. 
(2001) considered the case of non-spherical accretion from 
a hot accretion flow (ADAF) surrounding the neutron star. 
The hot ions heat the surface layer to produce a blackbody¬ 
like spectra with a high-energy Comptonized tail, highlighting 
the contribution of neutron star surface to the hard emission 
in LMXBs. Further investigations are required to determine 
how the surface emission is modified by the accretion flow. 
But in all cases, accretion onto the neutron star surface pro¬ 
duces a hardened blackbody-like spectrum and, at high ac¬ 
cretion accretion rates, a significant high-energy tail. 

We note that at luminosities above Lx ~ 10 36 erg s” 1 , cor¬ 
responding to groups 5-7, there is not a significant change in 
the hardness ratio or spectral parameters (photon indices, 
blackbody temperatures, etc.) for either the powerlaw or 
powerlaw plus blackbody models, despite over an order of 
magnitude change in luminosity. This ’plateau’ in spectral 
parameters could suggest that the emission, whether entirely 
nonthermal or a combination of thermal and nonthermal 
sources, above 10 36 erg s” 1 only changes in brightness with¬ 
out a substantial change in the accretion flow properties. 

4.3. Sub-Outburst Accretion and Crustal Heating 
Implications 

Accretion in the outburst state dominates the thermal evo¬ 
lution of transiently accreting neutron stars. During outburst, 
the accreted material compresses the upper layers of the neu¬ 
tron star, inducing pycnonuclear reactions in a process known 
as deep crustal heating (Brown et al. 1998). The heat gen¬ 
erated by the reactions is conducted throughout the neutron 
star and is partially radiated away at the surface, emerging 
as a thermal spectrum. From the time-averaged outburst ac¬ 
cretion rate, which is determined by the outburst history, one 
can estimate the neutron star’s quiescent luminosity and sur¬ 
face temperature due to deep crustal heating. Time-averaged 
accretion rates due to sub-outburst accretion, however, are 
less well-known as there are typically large uncertainties in 
the accretion events’ properties, such as duration, intensity 
and recurrence times as this type of low-level activity is diffi¬ 
cult to observe with all-sky X-ray monitors. 

The heating problem is particularly important in SAX 
J1750.8-2900 because it exhibits two forms of accretion out¬ 
side of outburst: sustained (> 1 month) accretion at 10 35 erg 
s' 1 and a low-level flare in quiescence, while also having one 
of the highest reported surface temperatures for a neutron 
star in quiescence (Lowell et al. 2012). The system warrants 
an investigation as to whether these forms of low-level accre¬ 
tion could be frequent enough to produce shallow heating in 
the crust and contribute to the high thermal luminosity. 

The intensity and type of accretion flow that manifests in 
quiescent transients has been a long-standing issue. A hard 
power law component often seen in quiescent neutron star 
spectra has commonly been attributed to ongoing accretion. 
It has been recently proposed that the hard, high-energy 
tail is produced by thermal bremsstrahlung emission from 
a radiatively-inefRcient accretion flow (Chakrabarty et al. 
2014). SAX J1750.8-2900, however, does not show signs of 
nonthermal emission in quiescence nor does it show signs 
of spectral variability, which is another sign of ongoing ac¬ 
cretion. The 2010 and 2012 quiescent luminosities were 
consistent (Lowell et al. 2012; Wijnands & Degenaar 2013) 
further supporting the claim that accretion is not contin¬ 
uous during quiescence in SAX J1750.8-2900, although the 
Wijnands & Degenaar (2013) luminosity has large uncertain¬ 
ties. 

While SAX J1750.8-2900 does not exhibit the character- 


isitic signs of continuous accretion in quiescence, the detec¬ 
tion of a flare (Wijnands & Degenaar 2013) suggests accre¬ 
tion is at least sporadic. Flaring activity of similar inten¬ 
sity and duration as in SAX J1750.8-2900 has been exhibited 
in several X-ray neutron star transient systems, including 
KS 1741-293 (Degenaar & Wijnands 2013) and XTE 1701- 
462 (Fridriksson et al. 2011). While we lack the statistics to 
monitor the progression of the flare in SAX J1750.8-2900, 
Fridriksson et al. (2011) studied the hardness-intensity dia¬ 
gram of a flare observed in XTE 1701-462 that reached a 
peak luminosity 20 times higher than the system’s quiescent 
level. Both the thermal and non-thermal fluxes increased, 
albeit differently, strongly suggesting the flare was accretion- 
powered. 

With our Swift XRT constraints on the accretion events 
(the quiescent flare and Lx = 10 34 — 10 36 erg s” 1 activ¬ 
ity) we can estimate their associated time-averaged accretion 
rates and compare to that of outburst, which is expected 
dominate the thermal evolution of the neutron star. We 
roughly approximated the flare as having a duration of 16 
days (the longest possible duration in order to calculate an 
accretion rate upper limit) with a constant flux (~4x 10” 12 
erg s” 1 cm -2 - our 0.5-10 keV unabsorbed flux estimate for 
the flare) and a (highly unconstrained) recurrence rate of 3 
times per year. Bolometric corrections are model and lumi¬ 
nosity dependent, but we use the ratio between the unab¬ 
sorbed bolometric and 0.5-10 keV quiescent fluxes reported 
by Lowell et al. (2012) as a correction estimate (B.C. = 1.2 
for the NSATM0S fluxes). After we convert the flare’s flux to 
a bolometric luminosity, we subtract the bolometric luminos¬ 
ity associated with the neutron star surface emission in qui¬ 
escence (i.e. the Lowell et al. (2012) quiescent luminosity). 
Assuming an accretion efficiency of e = 0.2, we find a time- 
averaged accretion rate of M = 1 x 10” 13 Mq per year (using 
LboI = eMc 2 ); even a recurrence rate of 10 flares per year 
only yields M = 3 x 10” 13 Mq per year which is nearly 2 
orders of magnitude smaller than our outburst averaged ac¬ 
cretion rates, M = 0.4 — 2.2 x 10” 10 Mg/year. 

We estimated the outburst M based on an outburst du¬ 
ration of 4 months with an occurrence rate of once every 
4 years. We have detected SAX J1750.8-2900 in outburst 
with luminosities between Lx = 0.4 — 2.5 x 10 37 erg s” 1 , 
corresponding to range of accretion rates previously stated. 
Lowell et al. (2012) computed an outburst mass accretion 
rate of 2 x 10” 10 Mq per year which is in agreement with 
our estimates. For SAX J1750.8-2900’s activity between Au¬ 
gust and October 2008, we estimated a maximum rate of 
M = 2 x 10” 12 Mq per year based on a duration of 60 days 
(August and September) with a typical Lx = 10 35 erg s” 1 
and a recurrence time of once per year. If this behavior is 
only associated with outbursts, then we’d expect a lower oc¬ 
currence rate (once per 4 years) and a lower mass accretion 
rate (~4x 10” 13 Mq). 

Although there are large uncertainties in our assumptions 
(recurrence times, durations, etc.) we find that sub-outburst 
accretion events lead to mass accretion rates at least an order 
of magnitude lower than rates associated with the outburst 
state. This supports the claim that the thermal evolution 
of the neutron star in SAX J1750.8-2900 is most heavily de¬ 
pendent on outburst accretion and that low-level accretion 
has little to no effect on the long term crustal temperature. 
Our conclusion is in agreement with Fridriksson et al. (2011), 
who performed a similar analysis based on flares exhibited by 
XTE J1701-462 and found that low-level accretion is unlikely 
to have a significant effect on the equilibrium surface temper¬ 
ature. We do not, however, address the immediate crustal 
heating effects due to small accretion events, such as flares, 
as they require much more detailed observations and calcula- 
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tions. 

5. CONCLUSIONS 

In our Swift XRT study of SAX J1750.8-2900’s X-ray 
spectral behavior over four years of activity, we have found 
that the source softens as its luminosity decreases between 
L x = 4 x 10 36 —»• 10 35 erg s _1 . This trend is consistent 
with the spectral softening observed in individual neutron 
star LMXB systems, as well as the cumulative softening be¬ 
havior exhibited by over one dozen transients between the 
luminosities of 10 36 and 10 34 erg s _1 . Our data is consistent 
with both the softening being due to a steepening power law 
towards lower luminosities and a thermal component becom¬ 
ing more apparent in the spectrum. The soft component may 
be associated with accretion onto the neutron star surface 
and has been definitively detected in other softening neutron 
star transients, which supports the interpretation that the 
thermal emission is driving the softening as the luminosity 
decreases. 
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